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Abstract 

At the first part of our communicate we show how specific umbral 
extensions of the StirHng numbers of the second kind result in new 
type of Dobinski-like formulas. In the second part among others one 
recovers how and why Morgan Ward solution of uncountable family of 
Tp- difference calculus nonhomogeneous equations A^/ — Lp'm the form 

extends to if)- Appell polynomials case automatically. Illustrative spec- 
ifications to g-calculus case and Fibonomial calculus case are made ex- 
plicit due to the usage of upside down notation for objects of Extended 
Finite Operator Calculus . 
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INTRODUCTION 



The upside down convenient notation for objects of Extended Finite Op- 
erator Calculus (EFOC) recently being developed and promoted by the 
present author ^ |21 El 13 is to be used throughout the whole exposition of 
both parts. 



The EFOC is an implementation of operator approach of Gian Carlo 
Rota to various by now formulations of extended of umbral difference calculi. 
For abundant references - thousands of them - see and ^ [21 El E] • 

In ^ 121 El S El and here also ip denotes a number or functions' sequence 
- sequence of functions of a parameter q. ip constitutes the Morgan Ward [6] 
extension (see: after then Viskov [7]) of (^)n>o sequence to quite arbitrary 
one (the so called - "admissible" in [8] and after then in [1,2, 3, 4, 5]). The 
specific choices are for example : Fibonomialy-extended sequence {^)n>o 
{{Fn) - Fibonacci sequence ) or just "the usual" ^-sequence (^)n>o or the 
famous Gauss g-extended {:;^)n>o admissible sequence of extended umbral 

operator calculus, where Ug = and Ugl = nq{n — 1)^!, Og! = 1 - see more 
below. With such type extension we frequently may in a "-;/'- mnemonic" way 
repeat reasoning and this what was done by Rota (see at first (FOC) of 
Rota in [9]). 

The simplicity of the first steps to be done while identifying general prop- 
erties of V'-extended objects consists in writing objects of these extensions 
in mnemonic convenient upside down notation ^ [21 E1I5 ) CHI which is 
here down introduced as follows: 

— ■ = n^,nA = n^{n-l)^\,n>Q,x^ = — , (1) 

x^ = x^{x-l)^{x-2)^...{x-k + l)^ (2) 
ihix — l)ibix — 2)...ib{x — k) 

If one writes the above in the form x^ = ^^^^^ = ^{x) = = , 
one sees that the name upside down notation is legitimate. 

You may consult ^ [21 [31 12 ISl and [TU] for further development and use 
of this notation . 
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In the first part of our communicate umbral extensions of the Stirhng 
numbers of the both kinds are considered and the resulting new type of 
Dobinski-hke formulas are discovered [10]. These extensions naturally en- 
compass the well known g-extensions . The fact that g-extended Stirling 
numbers giving rise to the umbral g-extended Dobinski formula interpreted 
as the average of powers of random variable Xq with the g-Poisson distribu- 
tion and are equivalent by re-scaling with the other Comtet -like g-extended 
Stirling numbers [10] - singles out the g-extensions itself which appear to 
be a kind of bifurcation point in the domain of umbral extensions . The 
further consecutive umbral extensions of Carlitz-Gould g-Stirling numbers 
are therefore realized in [10] in a two-fold way. 

In the second part of our communicate one displays jJJ how and why 
Morgan Ward solution |B] of ip- difference calculus nonhomogeneous equation 
A^/ = in the form 



recently proposed by the present author ^21" extends here now to tp- Appell 
polynomials case - almost automatically. Illustrative specifications to q- 
calculus case and Fibonomial calculus case El were already made 
explicit in ^2] exactly due to the of upside down notation for objects of the 



The First Part - Stirling numbers extensions and Dobinski- 
like formulae 

In this part of our communicate we follow [10] and we refer Reader for 
further details to consult [10]. At the start let us recall that the two standard 
^SIEIEI g-extensions Stirling numbers of the second kind might be defined 
as follows: 




EFOC . 



n 




(4) 



where Xq = and x^ = Xq{x — l)q...{x — k + l)q , which corresponds 
to the ip sequence choice in the (?-Gauss form (;;^)n>o 
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and (;~-Stirling numbers 

n 
fe=0 

where Xk{x) = x{x - lq){x - 2q)...{x -[k- l]q) 

For these two classical by now g-extensions of Stirling numbers of the 
second kind - the " (/-standard" recurrences hold respectively: 

rr}>i:(:)/-{.!j:-^o.*^i- 

From the above it follows immediately that corresponding g-extensions 
of Bn Bell numbers satisfy respective recurrences: 

Bq{n + l)=Y,{']] q'Bg{iy,n>0, 

1=0 ^'^9 
1=0 ^^^1 

where 

k=0 ^ 

Note that both definitions via (4) and (5) equations consequently cor- 
respond to different (/-counting jJB] . 

For applications to coherent state phisics see [17] and references therein. 

With any other choice out of countless choices of the tp sequence the equa- 
tion (5) becomes the definition of ip'^- Stirling (vide " Fibonomial-Stirling" ) 

numbers of the second kind|]^| and then ip'^-Bell numbers B!^{ip) are 

defined as sums as usual - where now Xk{x) in (5) is to be replaced by 
fpk{x) = x{x — — 2^)...(x — [k — l]ip). These ip^- Stirling numbers of 

the second kind recognized for q case properly as Comtet numbers in Wag- 
ner's terminology [10] satisfy familiar recursion and are given by familiar 
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formulas to be presented soon. 



The extension of definition (4) of the q- Stirling numbers of the second 
kind beyond this g-case i.e. beyond the ij} = (^)n>o choice is not that 
mnemonic and the " behavior " of the naturally expected recursion under 
extension comprises quite a surprise (see: appendix 2.2 in [10]). 

Therefore further consecutive umbral extension of Carlitz-Gould g-Stirling 
numbers ^-iid is realized two-fold way [10]. 

The first " easy way" consists in almost mnemonic sometimes replacement 
of q subscript by ip because we learn from 16., that the equation (5) defines 
as a matter of fact the specific case of Comtet numbers jl6| i.e. we define 
V'-extended Comtet-Stirling numbers of the second kind as follows: 

n 
A:=0 

where = x{x - l^){x - 2^)...{x - [k - 1]^). 

As a consequence we have "for granted" : 

where |q| = (5n,0) =0, k > n; and the easy derivable 

recurrence relations for ordinary generating function now read 

where 

n>0 

from where as a consequence we arrive to what follows: 



k 

Gk^{x)= T wi ^ wi y ^ ' 
'I' (1 - l^x)(l - 2^a;)...(l - fc^x) 
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and adapting reasoning from ^H] we derive the following explicit formula 



Expanding the right hand side of (9) results in another explicit formula 
for these ■i/'-case Comtet numbers i.e. we have 



{^} ^ X] (n)i/.(^2)v---(^n-fc)v; n,k>0. (11) 

l<ii<i2<---<in-k<k 

or equivalently (compare with [16]) 

{1)1= E -^k>0. (12) 

di+d2+---+rffe="~fci di>0 

With help of -(/'"-Stirling numbers of the second kind being defined 
equivalently by (6) , (7), (11) or (12) we define now ^~-Bell numbers in 
a standard way 

n 



k=0 



The recurrence for B!^{ip) is already quite involved and complicated 
for the g-extension case (see: the first section in [10])- and no acceptable 
readable form of recurrence for the ^-extension case is known to us. 

Nevertheless after adapting the the corresponding Wilf's reasoning from 
jlSj we derive for two variable ordinary generating function for | ^ | Stirling 

numbers of the second kind and the ^-exponential generating function for 
B'^{tp) Bell numbers the following formulae 

C;(x,y) = ^^:(V,y)x", (13) 

n>0 

where the ip- exponential-like polynomials A'^{'ilj,y) 



do satisfy the recurrence 



k=0 
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Ki'^^y) = + ^V']^n-i(V',y) « > i, 

hence 

A-{^,y) = [y{l + d^ri, n>0, 

where the Hnear operator acting on the algebra of formal power series is 
being called (see: [1,2,3,4,5] and references therein) the " '0-derivative" and 

The ■(/'-exponential generating function F^{x) = X^n>o -^n (V')^ 

for F~ Bell numbers - again - after cautious adaptation of the method 

from the Wilf's generatingfunctionology book [18] we get only a little bit 

involved formula 



r>0 

where (see: [6,7,1,2,3,4]) 



5^;(^) = E<V',r)^^ (14) 

r>0 



X 



n>0 

while 

°° ( — '\)k—r 

and for ^-extension the Dobinski like formula here now reads 

B-W=Y,e{Ar)% (16) 

In the case of Gauss g-extended choice of (;^)n>o admissible sequence 
of extended umbral operator calculus we have then 



^(-•)^E^.-® (17) 

k=r ^ 



and the new Q'~-Dobinski formula is given by 



r>0 
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which for = 1 becomes the Dobinski formula from 1887 [19]. 

As for the problem of how eventually one might interpret the '(/'~-Dobinski 
formulae (16) and (18) in the Rota-like way see: [10]. 

In [10] you find also details in support of our conviction that g-extensions 
seem appear as a kind of bifurcation point in the domain of umbral exten- 
sions. 



For the discussion of the other way to arrive eventually to another type 
of Dobinski formula - via an attempt to -(/'-extend the equation (4) - we again 
refer the Reader to [10] as this is quite a longer story with a surprise. 



The parallel treatment of the Comtet 



n 
k 



Stirling numbers of the 



first kind is now not difficult. 

Namely we define as in [10] the tp'^ Stirling numbers of the first kind as 
the coefficients in the following expansion 



r=0 



k 
r 



(19) 



where - recall ipk{x) = x{x — l^){x — 2^)...{x — [k — 1]^); Prom the 
above we then get 



i:[:i,{;}>^«- 



r=0 



(20) 



Another (expected Whithney numbers of the first kind) Stirling num- 
bers of the first kind [10] are defined as follows: 



r=0 



k 
r 



(21) 



where - now ip^{x) = x{x + l^)(x -|- 2^)...(x +[k — 1]^); More on that 
see [10]. 



The Second Part Appell polynomials' solutions of the 
Q{d^)- difference nonhomogeneous equation 

The second part is based on [11,12] to which we refer for more details. 
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At first recall the simple fact to be used in what follows. 

Proposition 0.1. Q{d^) is a ip- delta operator iff there exists invertible 
S € such that Q{d^) = d^S. 

Formally: "5 = Q/d^" or = In the sequel we use this 

abbreviation Q{d^) = Q. 

ijj- Appell or generalized Appell polynomials {^n(a^)}„>o defined ac- 
cording to 

d^An{x) = n^An-i{x) (22) 
and they naturally do satisfy the ip - Sheffer- Appell identity |H1 ^ 

A4x+^y) = j2('') ^s(y)x"-^ (23) 

tp- Appellor generalized Appell polynomials {An{x)}^^Q are equivalently char- 
acterized via their ip- exponential generating function 

Y.z^^^ = Aiz)exp^{xz}, (24) 

where A{z) is a formal series with constant term different from zero - here 
normalized to one. 

The ip- exponential function of V'- Appell- Ward numbers A^ = An{0) is 

Y,z^—,=A{z). (25) 

Naturally ^p- Appell {An{x)}^yQ satisfy the ip- difference equation 

QAn{x) = n^x""^; 'n > 0, (26) 

because QAn{x) = QS~^x"' = Q{d^/Q)x"- = d^x"^ = n^x"~^ ; n > 0. 
Therefore they play the same role in Q{d^)- difference calculus as Bernoulli 
polynomials do in standard difference calculus or ip- Bernoulli- Ward polyno- 
mials (see Theorem 16.1 in jj| and consult also ^2j) in ^-difference calculus 
due to the following: The central problem of the Q{d^) - difference calculus 
is: 

Q{d^)f = ^ =?, 
where f,(p- are for example formal series or polynomials. 
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The idea of finding solutions is the ip- Finite Operator Calculus mSlEHU 
standard. As one knows ^El) ^-ny V'" delta operator Q is of the form 

whereS G S^. Let Q{d^) = J2k>i I^.^ip^ 1^ 0. Consider then Q{d^ = 
d^S) with S = E,>o T^T^d'^ ^ tJo ^9'; so = qi^ 0. 

We thus have for = A - call it: V" Appell operator - the obvious 
expression 

Now multiply the equation Q{d^)f = hy A = ^n>o n^.^^ thus 
getting 

^^/ = E^^^"^' V'("^ = 5^^^"-'^- (27) 
The solution then reads: 

f{x) = Y,—^^^''-^\x)+ [ vix)+pix), (28) 

where p is "Q(5^)- periodic" i.e. Q{d^)p = 0. Compare with 
for "+^1- periodic" i.e. p(x +^ 1) = p{x) i.e. A^p = 0. Here the 
relevant ^ip - integration V^l^^) is defined as in [T. We recall it in brief. Let 
us introduce the following representation for " difference- ization" 

= h^,dQ ; n^x"-"-^ = n^x""-^; n > 1, 

where Oqx^ = i.e. do is the q = Jackson derivative, do is identical 
with divided difference operator. Then we define the linear mapping 
accordingly: 



X- = I x-i- I x" = - — ^— x"+^ n > 



n^y (n + 1)^ 



where of course 5^ o L = id- 
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1 Examples 

(a) The case of ip- Bernoulli- VFard polynomials and A^- difference calculus 
was considered in detail in J2] following 6 . 

(b) Specification of (a) to the Gauss and Heine originating g'-umbral cal- 
culus case El in El m was already presented in 12[. 

(c) Specification of (a) to the Lucas originating FFOC - case was also pre- 
sented in J2] (here: FFOC=Fibonomial Finite Operator Calculus), 
see example 2.1 in P]). Recall: the Fibonomial coefficients -already 
known in 19-th century to Lucas [20]- are defined "binomially" as 



where (-F„- Fibonacci numbers in up-side down notation: np = F^ ^ 0, 
up^- = np{n — l)p{n — 2)p{n — 3)f ■ ■ - ^plp; Opl = 1; 



We shall call the corresponding linear difference operator dp; dpx"^ = 
npx"'~^; n > the F-derivative. Then in conformity with ^ and 
with notation as in [T]- jl3l HH] one has: 



for the corresponding [3,1,4] generalized translation operator E°'{dp). 
The V'- integration for the moment is still not explored F- integration 
and we arrive at the F- Bernoulli polynomials unknown till now. 

Note: recently a combinatorial interpretation of Fibonomial coeffi- 
cient has been found Jl] by the present author. 

(d) The other examples of Q{d^)- difference calculus - expected naturally 
to be of primary importance in applications are provided by the pos- 
sible use of such ijj- Appell polynomials as: 

• ■i/'-Hermite polynomials {Hn^ip}^-^^- 
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k>0 



k,i,l 



n > 0; 



tp - Laguerre polynomials {in,i/;}„>o 



n 



n / ^ 

k=l ^ ^ 

n 



For q = 1 m (/-extended case one recovers the known formula : 



Ln,,=l{x) = E (-1)* 
k=l 



kn'f n-1 



^=9! \ k-l 



X . 



References 



[1] A.K.Kwasniewski Main theorems of extended finite operator calculus 
Integral Transforms and Special Functions, 14 No 6 (2003): 499-516 

[2] A.K.Kwasniewski On Simple Characterizations of Sheffer psi- 
polynomials and Related Propositions of the Calculus of Sequences , Bul- 
letin de la Soc. des Sciences et de Lettres de Lodz, 52 Ser. Rech. Deform. 
36 (2002):45-65, ArXiv: |math.CQ /0312397| . 

[3] A.K.Kwasniewski Towards ip-extension of Finite Operator Calculus of 
Rota Rep. Math. Phys. 48 No3 (2001):305-342 

[4] A.K.Kwasniewski On extended finite operator calculus of Rota and 
quantum groups Integral Transforms and Special Functions 2 No4 
(2001):333-340 

[5] A.K.Kwasniewski First contact remarks on umbra difference calculus 
references streams Bulletin de la Soc. des Sciences et de Lettres de 
Lodz, in press; ArXiv: math.CO/0403139 vl 8 March 2004 

[6] M. Ward: A calculus of sequences, Amer.J.Math. Vol.58, 1936, pp. 255- 
266 



12 



O.V. Viskov: Operator characterization of generalized Appel polynomi- 
als Soviet Math. Dokl. 16, 1521 (1975). 

G. Markowsky, Differential operator and theory of binomial enumera- 
tion, Math. Anal. Appl. 63 (1978), 145-155. 

G. -C. Rota: Finite Operator Calculus, Academic Press, New York 1975. 

A.K.Kwasniewski On umhral extensions of Stirling num- 
bers and Dobinski-like formulas submitted for publication 
|ArXiv:ma th.CO/0411002 30 October 2004 

A.K.Kwasniewski ip-Appell polynomials' solutions of an umbral differ- 
ence nonhomogeneous equation Bulletin de la Societe des Sciences et 
des Lettres de Lodz (54) Serie: Recherches sur les Deformations 45 
(2004), 11-15. ArXiv: math.CO/0405578 30 May 2004 

A. K. Kwasniewski: On basic Bernoulli- Ward polynomials , Bulletin de 
la Societe des Sciences et des Lettres de Lodz (54) Serie: Recherches 
sur les Deformations 45 (2004), 5-10. |ArXiv:math.CO/0405577| 30 May 
2004 

A.K.Kwasniewski Fibonomial cumulative connection constants Bulletin 
of the ICA 44 (2005), 81-92. ArXiv: math.CO/0406006 1 Jun 2004 

A.K.Kwasniewski: More on combinatorial interpretation of the fibono- 
mial coefficients Bulletin de la Soc. des Sciences et de Lettres de Lodz 
(54) Serie: Recherches sur les Deformations 44 (2004), 23-38. ArXiv: 
|math.CO/0402344 vl 22 Feb 2004 

S.C. Milne A q-analog of restricted growth functions, Dobinski's equal- 
ity, and Charlier polynomials , Trans. Amer. Math. Soc. 245 (1978) 
89-118 

Carl G. Wagner Partition Statistics and q-Bell Numbers (q = -1), Jour- 
nal of Integer Sequences 7(2004) Article 04.1.1 

J.Katriel, M. Kibler Normal ordering for deformed boson operators and 
operator-valued deformed Stirling numbers J. Phys. A: Math. Gen. 25 
(1992): 2683-26-91 

H. S.Wilf Generating functionology Boston: Academic Press, 1990 

G. Dobinski Summierung der Reihe S .... fiir m = 1, 2, 3, 4, 5, .... 
Grunert Archiv (Arch. Math. Phys) 61, 333-336, (1877) 

[20] E. Lucas, Theorie des fonctions numeriques simplement periodiques, 
American Journal of Mathematics 1 (1878): 184-240 



13 



